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We consider the interpretation of a pair of kHz Quasi Pe- 
riodic Oscillations (QPOs) in the Fourier spectra of two Low 
Mass X-Ray Binaries, Sco X-1 and 4U1608-52, hosting an old 
accreting neutron star. The observed frequency difference of 
these QPOs decreaseas as their frequency increases, contrary 
to simple beat frequency models, which predict a constant 
frequency difference. We show that the behaviour of these 
QPOs is instead well matched in terms of the fundamental 
frequencies (in the radial and azimuthal directions) for test 
particle motion in the gravitational field of the neutron star, 
for reasonable star masses, and nearly independent of the star 
spin. The radial frequency must be much smaller than the az- 
imuthal one, testifying that kHz QPOs are produced close to 
the innermost stable orbit. These results are not reproduced 
through the post-Newtonian (PN) approximation of General 
Relativity (GR). kHz QPOs from X-ray binaries likely provide 
an accurate laboratory for strong field GR. 

PACS numbers: 97.60.Jd, Q4.80.Cc, 97.8Q.Jp 



Introduction. The kHz QPOs in the X-ray flux of 
some 20 old accreting neutron stars in Low Mass X- 
Ray Binaries (LMXRBs) are the first involving timescales 
close to the dynamical ones in the vicinity of a collapsed 
star The frequency of these QPOs (from ^ 0.33 to 
~ 1.23 kHz) displays changes which are often correlated 
with the instantaneous accretion rate. A common phe- 
nomenon is the presence of a pair of kHz QPOs (with 
centroid frequencies vi and i'2) drifting towards higher 
frequencies while keeping their difference Av = 1/2 ^ 
approximately constant. The range of Av is also quite 
narrow across different sources (A;^ k, 250 — 360 Hz). 
During type I bursts from six of these sources, a nearly 
coherent signal at a frequency of k, 330 — 590 Hz was also 
detected. In a few cases the latter frequency is consistent, 
to within the errors, with the frequency separation Av, 
or twice its value 2Az/. 

These results have been interpreted in terms of beat 
frequency models, BFMs (magnetospheric and sonic 
point) 1^ where inhomogeneities at the inner disk bound- 
ary are accreted at the beat frequency between the lo- 
cal Keplerian frequency, vki and the neutron star spin 
frequency, Vspin- i^K increases in response to increas- 
ing accretion rates and so does the QPO beat frequency 
i^b = vk — Vspin- In these BFMs, the lower frequency kHz 



QPOs are interpreted as the beat frequency (y\ — vt) 
and the higher frequency kHz QPOs as the Keplerian fre- 
quency {1/2 = vk) at the inner disk boundary. Therefore, 
Av yields the neutron star spin frequency [vspin — Av). 
A signal at a frequency close to Av (or 2Az^) is occasion- 
ally detected, but only during type I bursts. The inferred 
i^spin spans a remarkably small range (^^ 250 — 360 Hz), 
far from the breakup limit of any neutron star model. 
A direct consequence of the BFM interpretation, such a 
clustering of spin frequencies across LMXRBs with very 
different average mass transfer rates, orbital parameters, 
evolutionary histories and magnetic field strenghts is sur- 
prising and calls for a separate explanation 

Recent findings, however, undermine the applicabil- 
ity of simple BFMs to a growing number of LMXRBs 
with kHz QPOs. First, the separation Ai/ decreases by 
~ 80 - 100 Hz as V2 increases from ~ 800 to ~ 1100 Hz 
in Sco X-1 |] and 4U1608-52 g. Owing to poor statis- 
tics, a similar variation of Ai/ in other sources would 
have remained undetected the behaviour of Sco X-1 
and 4U1608-52 might therefore be the rule. Second, in 
4U1636-536 the frequency of the signal during bursts is 
inconsistent with Av or its second harmonics Q. 

We propose here a substantially different scenario in 
which both kHz QPOs originate from the motion of mat- 
ter inhomogeneities in the vicinity of the neutron star: 
as in BFMs, the higher frequency QPOs at V2 are pro- 
duced by the ^-motion (i.e. the Keplerian motion) of 
inhomogeneities orbiting the inner disk boundary, while 
the lower frequency QPO signal at vi originates from the 
periastron precession at the inner edge of the accretion 
disk. Strong-field effects play a major role in determin- 
ing the latter frequency. If this interpretation were con- 
firmed, kHz QPOs would provide an unprecedented test 
for strong-field general relativity, GR. 

The model. We assume that the motion of the disk 
matter is dictated by the star's gravity alone. The gravi- 
tational field close to a neutron star is dominated by the 
Schwarzschild terms of the metric, with small corrections 
due to the neutron star rotation. For nearly circular or- 
bits (eccentricity e ~ 0) the coordinate frequency in the 
radial direction vs,r as measured from an observer at in- 
finity is given by 



i's,r = J^s,4l-6MA)i/2 



(1) 
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where r is the coordinate distance and M the mass of 
the compact object (throughout this paper we use geo- 
metrical units G — c — 1). The (/)-frequency retains its 
Keplerian value vs.^ = vk = {M j A.-K^r')^!'^ . 

Eq. (1) shows that approaching the marginally sta- 
ble orbit (r„is = 6M), vg ^ vanishes, so that the pre- 
cession frequency of the orbit's periastron fper = — 
Vr = 1^5,0 — vs -r approaches the (/)-frequency. We pro- 
pose here that, while the QPOs at V2 arise from the 0- 
frequency of matter orbiting close to the neutron star 
(as in BFMs), the QPOs at v\ result from the strong- 
field periastron precession of the same matter. Therefore 
b.v = vi - vx = vs.tf, - {vs.,4> - vs,r) = vs.r- The curves 
in Fig.l show vs,r as a function of i/s,</) for selected val- 
ues of the neutron star mass, the only free parameter in 
the equations above. It is easy to see that vs^r decreases 
after a maximum of i's.r{'>nax) — 358(2Mq/M) Hz for a 



-frequency of vs.ip — 2vs^r{'max) ~ 0.65 I's^tpi'f 



The 



measured Av vs. 1^2 for ten LMXRBs is also plotted. It is 
apparent that for neutron star masses in the 2 Mq range, 
the simple model outlined above is in qualitative agree- 
ment with the measured values, including the decrease of 
Ai' for increasing 1^2 seen in Sco X-1 and 4U1608-52. 

Developing a quantitative model involves allowing for 
finite orbital eccentricity and neutron star rotation. In 
order to calculate the 0- and r-frequencies of bound orbits 
in the Schwarzschild metric, we define /i = M/a{l — e^), 
where a is the elliptic orbit semi-major axis. The radial 
frequency is given by S 



where 

/(e,M) = 



(l_e2)3/2((2^-l)2_4^2g2)l/2 ' 

cfe[l -2/x(3-f ecosa;)]-i/2 
(1 -I- ecosx)2[l — 2^(1 -I- ecosa;)] 



(2) 



(3) 



and Vk = {M/iTT^a^Y^^ . The azimuthal frequency 1^5 ^ 
is defined by an identical formula, except that the inte- 
gration in X extends to Xc, implicitly defined by 



27r-\/l — 6/i + 2/ie 



(l-Fsin^ 7)1/2 



(4) 



with k = 4/ie/(l — 6/i + 2/ie). 

The effects of the neutron star rotation are twofold. 
The gravitomagnetic effects are estimated here at the 
first order of the angular momentum in the Kerr met- 
ric, i.e. the Lense-Thirring order. We have |^ ^lt,cI} = 
-VLT,nodl^ and VLT.r = 'ivLT.nodl^ , whcrc 



VLT,nod = 47raV^(l - e 



(5) 



is the Lense-Thirring nodal precession frequency, a! — 
2'Ki'spinI /M the specific angular momentum and / the 
moment of inertia of the rotating star. 



The neutron star rotation-induced oblateness gives rise 
to a quadrupole component of the gravitational field. 
For orbits close to the equatorial plane, the quadrupole 
corrections are VQ,,f, = —VQ,nod{^ — e^Y^"^ and vq^^ = 
i^Q,nod[2 - (1 - e'^y/^] where 



'^Q,nod = -3z/K*2(2Afr2)-i(l - e^)- 



(6) 



is the quadrupole nodal precession frequency and $2 the 
star's quadrupole moment. Detailed numerical calcula- 
tions and approximate formulae for VLT.nod and VQ^nod 
are given in 1 10| for a variety of equations of state, EsOS. 
The corresponding corrections to the (j), r and perias- 
tron frequencies are small (up to tens of Hz) also for 
fast spinning neutron stars. Note that nodal precession 
arises exclusively from the relevant Lense-Thirring and 
quadrupole terms. 

The (j)- and r-frequencies obtained by summing the 
components discussed above were fit to the 1/2 and Av 
frequencies measured in Sco X-1. Initially M and e were 
treated as parameters of the fit, while different frequen- 
cies were calculated within each model by stepping a, 
with Vspin = 0. These fits had e ~ and were similar to 
those in Fig. 1. Allowing for neutron star rotation (in- 
cluding the Lense-Thirring and quadrupole terms) gave 
improved fits only for very stiff EsOS, such as L (c/. |l^] 
and references therein), and large counterrotating spin 
frequencies. For example a model consisting of nearly 
circular orbits with Vspin = —520 Hz, and a 1.8 Mq neu- 
tron star with EOS L gave a fit similar to that shown in 
Fig. 2. While viable, this model is not satisfactory, be- 
cause, if the neutron stars in LMXRBs have been spun 
up by accretion torques, one expects the disk matter to 
rotate in the same direction as the neutron star. 

Since there is no a priori reason why matter in the in- 
nermost disk region should maintain a constant e as a 
varies, we stepped e in order to obtain different frequen- 
cies. The neutron star equator, marginally stable orbit 
or magnetospheric radius (depending on which is largest) 
provide a natural inner boundary to the accretion disk 
and the periastron of its innermost orbits. With this 
geometry in mind we introduced the periastron distance 
Tp = a(l — e) as the second parameter in the fit. The best 
fit obtained in this way for a non-rotating neutron star 
is shown in Fig. 2. The agreement with the measured 
frequencies is remarkable. The variation of e with is 
also shown. The model has, for the two free parameters. 



M ~ 1.90 M(T, and 



6.25 M, a value close to the 



marginally stable orbit. 

When Vspin is allowed a finite value, slightly improved 
fits are obtained. For acceptable values of the neutron 
star spin frequency, however, these would be hard to dis- 
tiguish from the best fit in Fig. 2. For a (corotating) 
Vspin of 300 and 600 Hz and EOS L, the best fit para- 
maters would be M ~ 1.94 and 1.98 M©, rp ~ 6.18 
and 6.17 M, respectively. This shows that the effects 
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induced by the neutron star rotation through the Lense- 
Thirring and quadrupole terms are quite smaU, though 
non-neghgible; they cannot be currently disentagled from 
the Schwarzschild terms without a serious risk of overin- 
terpreting the data and model. While consistent with 
it, our model does not provide independent evidence that 
Sco X-1, and likewise other LMXRBs, contain a fast spin- 
ning neutron star. 

We have shown that the model we propose accurately 
reproduces the observations with a minimum number of 
free paramaters (two, M and r^) the values of which are 
well within the expected range. The fits described above 
(plus a large number of others that we tried), together 
with analytical estimates, show that for no acceptable 
value of the neutron star parameters (mass, spin, EOS) 
the Lense-Thirring and quadrupole corrections alone can 
produce a large enough decrease of Vr over the relevant 
range of i'^. Also, one may wonder whether the ac- 
curacy of our fits reflects the properties of the strong 
field Schwarzschild metric, or some lower order expan- 
sion would reproduce the observed frequencies compara- 
bly well. So we carried out the calculations as in p| with 
the approximate PN metric for a spherical star |ll ] , and 
repeated the fits described above. The modelling of the 
Sco X-1 data was very poor, the best fits lying outside 
the range of Fig. 2. We conclude that the matching of 
the model and the data is to be ascribed mainly to the 
strong-field Schwarzschild terms. 

A simple argument shows also that the strong field 
regime applies. In the interpretation we propose, i^r, 
often referred to as the epicyclic frequency, is <C v,p. 
This ordering is unusual, the epicyclic frequency in most 
cases being the larger. However, if one remembers that 
47r^^'^ = d'^Vef f /dr"^ , where is the effective poten- 
tial, one sees that the above ordering is natural for or- 
bits close to the marginally stable one, for which the 
effective potential notoriously has an inflection point, 
d^Veff/dr^ — 0. Thus the relevant orbits must be close 
to the innermost stable orbit. 

Discussion. The higher frequency kHz QPOs at 1^2 
are interpreted in terms of the orbital (/)-frequency of mat- 
ter at the inner edge of the accretion disk both within 
BFMs and the model presented here. The similarity in 
the highest value of 1/2 1-1 — 1-2 kHz) across different 
sources and, in a few cases, the saturation of 1/2 at its 
highest value despite an increasing source flux have been 
interpreted by several authors as signs that the inner disk 
edge reaches the marginally stable orbit . 

If, as predicted in our model, the lower frequency kHz 
QPOs arise from the periastron precession of the same 
matter that produces the higher frequency kHz QPOs, 
then the exciting perspective opens up of probing directly 
the predictions of GR in the vicinity of the marginally 
stable orbit. Moreover, results which are difficult to rec- 
oncile with BFMs are readily interpreted. First, the de- 
creasing Ai/ for increasing 1^2 observed in Sco X-1 0] and 



4U1608-52 (and compatible with measurements in all 
other kHz QPO LMXRBs) results naturally from the ex- 
pected strong field behaviour of t'r in the vicinity of the 
marginally stable orbit. Second, the fact that in all kHz 
QPO sources Af is around 250 — 360 Hz is naturally in- 
terpreted in terms of the characteristic broad maximum 
of Vr versus Vf/, apparent in Fig. 1 provided the neutron 
star is 2 M0. The predicted Vr and depend only 
weakly on Vspin- Therefore, the measurement of Vspin in 
QPO LMXRBs has to rely only on coherent pulsations 
(so far detected only in ISAX 1808.4-3658 Q) or, per- 
haps, the signals observed during type I bursts from (low 
luminosity) Atoll sources [Q . The range of neutron star 
spins (~ 330 — 590 Hz) inferred from Uburst is wider and 
shifted towards higher values, compared to that obtained 
within BFMs from Av. This helps reconciling the spin 
period distribution of LMXRBs with that of millisecond 
radio pulsars. Cyg X-2 and GX 17+2, somewhat pe- 
culiar type I bursters, might soon provide the first spin 
measurements in (high luminosity) Z-sources. 

The model described here requires a small eccentric- 
ity of the blobs' orbits at the inner disk edge in order 
to match in detail the Sco X-1 kHz QPOs, without in- 
voking large counterrotating neutron star spins. Such an 
eccentricity might originate from a variety of effects such 
as the amplification of the blobs' radial oscillations due 
to a resonant interaction with the rotating magnetic field 
of the neutron star [|l^ or quasi-elastic impacts of blobs 
with the neutron star magnetospheric boundary or the 
boundary layer at the neutron star surface. If eccentric- 
ity variations drive the observed changes of vi and 1^2, 
then the correlation of vi and 1^2 with the mass accretion 
rate inferred from observations requires an anticorrela- 
tion of the eccentricity and M ; the mechanism for this is 
presently unknown, but one may speculate that in low- 
eccentricity encounters of the blobs with the edge of the 
precipice there is more time for matter to fall over. 

A modulation at the periastron precession frequency 
can be produced through the usual mechanisms, namely 
orbit-modulated emission by self-luminous blobs or par- 
tially occulting blobs at the innermost disk boundary 
(c/.also fl^). In the first case Doppler beaming and 
blueshift enhance emission by the blobs when the perias- 
tron velocity (the highest along the orbit) is closest to the 
line of sight; in the second the fraction of radiation from 
close to the neutron star that is occulted is a maximum 
whenever occultation is produced by a blob at periastron, 
where a blob subtends the largest solid angle relative to 
the neutron star. Both situations repeat themselves with 
frequency — Vr, the periastron precession frequency. 
The harmonic content of the signal(s) observed at infin- 
ity might be affected by strong gravitational lensing in 
the vicinity of the neutron star |l^ . 

Finally we note that the present model has also impli- 
cations for the recent suggestion that the ~ 20 — 35 Hz 
oscillations observed in several kHz QPO sources of the 
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Atoll group, and perhaps also the ^ 15 — 60 Hz hori- 
zonthal branch oscillations of Z-sources, originate in the 
(Lense-Thirring dominated) nodal precession of blobs 
lifted off the equatorial plane at the inner disk boundary 
| p6[ . By using the spin frequencies inferred from Vburst 
(as opposed to Av in BFMs) the predicted nodal preces- 
sion frequencies for Atoll sources are up to a factor of ~ 2 
higher and match the observations more accurately . 
On the contrary, in the application of the nodal preces- 
sion model to Z-type LMXRBs the (as yet unmeasured) 
spin frequency is to be regarded as a free parameter. 

In summary, the QPO phenomenon in LMXRBs is 
likely to reflect the fundamental frequencies of the rel- 
ativistic motion of matter in the vicinity of the neu- 
tron star. As discussed in the present Letter, if the 
^-frequency and pcriastron precession frequency are di- 
rectly observed as kHz QPOs, these yield an unprece- 
dented diagnostic of GR in the strong-field regime. 



FIG. 1. kHz QPO frequency difference Av versus higher 
QPO frequency 1/2 for ten LMXRBs (of. [4, 5, 6] and references 
therein). Error bars are not plotted for the sake of clarity. 
The models give the r- and 0-frequencies of matter in nearly 
circular orbit around a non-rotating neutron star, of mass 2.2, 
2.0 and 1.8 Mq. 

FIG. 2. kHz QPO frequency difference Av versus higher 
QPO frequency 1/2 in Sco X-1. The best fit model corresponds 
to the r- and <j!>-frequencies of matter orbiting a non-rotating 
1.90 M0 neutron star at a periastron distance of 6.25 M (17.5 
km). The line marked with e gives the orbital eccentricity 
(xlOOO) as a function of the ^-frequency. Note that e ap- 
proaches for the highest observed 1^2', however a similar fit 
would be obtained by forcing rp to a slightly smaller value 
such that e remains finite, though small, for any observed 1/2- 
The upper curve is obtained by using the same parameters as 
above for a neutron star spinning at 600 Hz and EOS L. 
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